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The evolution of large-scale density perturbations is studied in a stably stratified, two-dimensional 
flow governed by the Boussinesq equations. As is known, intially smooth density (or temperature) 
profiles develop into fronts in the very early stages of evolution. This results in a frontally 
dominated potential energy spectrum. The fronts, initially characterized by a relatively simple 
geometry, spontaneously develop into severely distorted sheets that possess structure at very fine 
scales, and thus there is a transfer of energy from large to small scales. It is shown here that this 
process culminates in the establishment of a k~ 5 ^ 3 kinetic energy spectrum, although its scaling 
extends over a shorter range as compared to the scaling of the potential energy spectrum. The 
establishment of the kinetic energy scaling signals the onset of enstrophy decay which proceeds 
in a mildly modulated exponential manner and possesses a novel self-similarity. Specifically, 
the self-similarity is seen in the time invariant nature of the probability density function (PDF) 
associated with the normalized vorticity field. Given the rapid decay of energy at this stage, the 
spectral scaling is transient and fades with the emergence of a smooth, large-scale, very slowly 
decaying, (almost) vertically sheared horizontal mode with most of its energy in the potential 
component — i.e. the Pearson-Linden regime. 



PACS numbers: PACS number 47.52.+j 
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I. INTRODUCTION 

Large-scale geophysical flows usually evolve under the constraints of stable stratification and rotation. Indeed, it 
is known that both of these constraints, individually and in concert, profoundly affect the motion of a fluid Q|. Here 
we restrict our attention to the stratified problem (see 0] for a review). In particular, we study the adjustment of a 
two-dimensional stably stratified Boussinesq fluid to imposed large-scale density perturbations. 

Starting with the work of Riley, Metcalfe & Weissman Q , there have been numerous studies of three-dimensional 
stably stratified flows in both forced and decaying scenarios that indicate a spontaneous generation of layered 
structures from initially isotropic fields (see for example 0,0,0). Indeed, the small Froude number limit of the 
three-dimensional Boussinesq equations (especially in the decaying case) indicates layerwise motion with weak 
vertical correlation Q , > 0L_ Further, amongst the variety of wave- vortex and wave- wave interactions possible in 
the three-dimensional case it appears as though the near resonant transfer into the slow (horizontal) wave 

modes via fast waves is responsible for the generation of these layered structures i. In fact the partitioning of 
energy between the fast and slow components in the general rotating stratified system is an outstanding problem in 
atmospheric dynamics (for example Ford et. al and the references therein). In light of this, it may be instructive 
to understand how energy re-distribution proceeds in the relatively simplified, two-dimensional system that no longer 
possesses a distinct vortex mode and presents a situation wherein one can study the interaction of wave modes in 
isolation. A related issue is the fact that potential vorticity is identically zero in two dimensions and a future goal 
is to understand if and how the constraint of potential vorticity conservation in three-dimensional flows alters wave 
interactions. 

With regards to the two-dimensional problem, detailed numerical work was carried out by Bouruet-Aubertot, 
Sommeria & Staquet |lC|.|ll| wherein a bounded domain with no-flow conditions was considered. Focussing on the 
instabilities suffered by prescribed flows, they showed that the wave-wave interactions result in a transfer of energy 
from large to small scales (see also Orlanski & Cerasoli 01 )■ Their analysis of imposed standing gravity waves 
indicated that these waves necessarily break after a finite time (depending on the initial amplitude) and the spectral 
redistribution of energy proceeds via the so-called wave-turbulence paradigm |13| . Here we lift the no flow condition 
by considering the problem in a periodic domain and pursue the time evolution, from a state of rest, of a smooth 
large scale density (or temperature) perturbation. As one might anticipate, the aforementioned wave-breaking is one 
stage in the overall scheme of things. 

Apart from geophysical motivatio ns, the system under consideration is part of an extended family of flows — the 
so-called dynamically active scalars p^|.|l5j|. In particular, statistical properties of the neutral and unstably stratified 
Boussinesq systems have been the subject of recent investigations E3;E3iEl an d the present work, wherein the 
system additionally supports waves, can be viewed as a continuation of these efforts. Further, given the possibility of 
a finite-time breakdown of regularity in a two-dimensional setting, these problems are of considerable mathematical 
interest and the reader inclined to pursue such matters is referred to Cordoba & Fefferman (and the references 
therein) for a recent overview. 

The remainder of the paper is organized as follows : we first provide the basic equations and conservation laws 
along with the setup of the numerical experiment. Next we touch briefly upon the initial stages of evolution that 
have been well documented in existing literature. We then proceed to the focus of the paper showing a k~ x scaling 
for the potential energy spectrum followed by the establishment of a transient k~ 5 ^ 3 scaling for the kinetic energy, 
the onset of nearly exponential decay of enstrophy and the establishment of an invariant normalized vorticity PDF. 
Due to the rapid decay of energy at this stage of the problem, the spectra are seen to gradually shift downwards, 
their scaling becomes less distinct, and finally we see the emergence of a very slowly decaying, vertically sheared, 
large-scale, predominantly horizontal mode — most of whose energy is trapped in the potential component — as 
elucidated by Pearson & Linden [2(j. We conclude by summarizing the various stages, discuss the fate of stratified 
adjustment in the inviscid limit and point out avenues for future research. 



II. THE GOVERNING EQUATIONS 

The equations governing two-dimensional stratified flow under the Boussinesq approximation are 
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In the above p = pb + p' ,p = Pb + p 1 and T = Tb + 9, where p& = po + In particular, /3(.z) = — £?z and the basic 
state is in hydrostatic balance, i.e. dpb/dz = —pb9- The temperature and density are related via p = p a [l + a(To — T)], 
and thus the basic state temperature profile is Tb = T + T(z) with T(z) = —p(z)/(ap ). System Q results from 
the assumption p' ,p{z) « p . As p Q1 X,g,a are constants we consider the perturbed fields (u,p' , p' ,9) in a periodic 
domain |38| . In vorticity-stream form Q reads 

Duo _ 89 D0 _ dtp 
~Dt ~ ~ a9 dx~ ' Dt ~ dx 

UJ = -Vy J;u = --; W = -. (2) 
Defining E = J D [u 2 + w 2 + ^j-6 2 ] , we see that (QJ conserves — which we refer to as the total energy of the system [3^ . 
By linearizing J5J it can be seen that the above supports gravity waves obeying the dispersion relation 



a(k) = ±V(Xgaj^ ; K 2 = k 2 x + k 2 z . (3) 

Note that the presence of a uniform gravitational field in (J2J breaks symmetry with respect to mirror images about 
the x-axis. Apart from this lack of complete parity invariance, in contrast to the usual two-dimensional scenario 
the vorticity in is not a Lagrangian invariant. This indicates that the familiar two-dimensional inverse cascade 
of kinetic energy is unlikely to be realised as it relies, amongst other thin gs, on the dual conservation of energy and 
enstrophy . Also, in comparison to the free-convective case (i.e. A = 0) [0|,[lg], functions of the form f(9) are not 
conserved by @ ■ Hence an immediate direct cascade of 9 2 is unlikely. Introducing dissipation via the usual kinematic 
viscosity and diffusivity (both taken to be 7) — i.e. 7V 2 w and jV 2 9 in — we have 
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So for 7 > 0, ^ < and |0J imply that all fields have to die out in the limit t — * 00. Of course, even though E is a 
monotonically decreasing function of time, the enstrophy (and higher moments of u>) can show a significant increase 
before finally decaying away. 



III. THE NUMERICAL EXPERIMENT 



To study the process of adjustment, we start with an initial state at rest and in hydrostatic balance to which is 
imposed a large scale 9 perturbation. Various mathematical forms for the initial condition were used to check the 
qualitative similarity in the response of the system, in particular the results presented correspond to to a Gaussian 
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bump (6(x, z, 0) = 0.5 exp{— 4[(x — tt) 2 + (z — 7r) 2 ]}) in the middle of the domain. System (|2J) is solved at a resolution 
of 750 x 750 in a 27r periodic domain using a pseudo-spectral method de-aliased by the two-thirds rule and a fourth 
order Runge Kutta time stepping scheme. Defining the Froude (Fr) and Reynolds (Re) numbers as Fr = u/NL 
and Re = uL/^ respectively, taking N to be an 0(1) quantity consider an initial 9 perturbation such that the total 
energy is also an 0(1) entity. Apriori we expect the following stages of evolution : as the fluid starts from rest, 
initially both Fr, Re < 1 — i.e. we have a laminar flow that is strongly constrained by stratification. Since 7 is 
small, as soon as a certain fraction of the potential energy is converted to kinetic energy implying max(u) ~ 1, we 
have Fr ~ 1, Re 3> 1 — i.e. a turbulent flow that is weakly constrained by stratification. Finally, because 7 > the 
fields must eventually decay so that we enter a diffusion and stratification dominated regime. In accordance with 
this scenario we set the parameters in the present simulation tobep o = l,-B = 0.1,a = 0.1,7 = 5xl0 -4 and TV = 1. 



A. The formation of fronts and their subsequent distortion 



Fig. 1 shows the evolution of energy — total, potential (J D ^d 2 ) and kinetic (J D [u 2 + w 2 ]) — and enstrophy with 
time. Focussing on the early stages of development shown in the first column of Fig. 2, as predicted from linear 
theory we see the generation of gravity waves that mediate the exchange of energy between the potential and kinetic 
components. Indeed, simulations with differing initial conditions and variations in the strength of stratification follow 
qualitatively similar paths though the quantitative partition of energy between the potential and kinetic components 
is not identical in all situations. Further, at very early times (t < 1.5s), all gradient fields are quite mild and the 
evolution is fairly inviscid. Along with this laminar evolution, Fig. 3 shows the emergence of frontal structures 
perpendicular to the vertical direction. Moreover, this frontal development is accompanied by increasing vertical 
shear and, even though this is not a steady flow, the accompanying decrease in Richardson number hints at the onset 
of instability [2j|. In fact, as noted in previous studies such a situation necessarily leads to wave- breaking ^(j. In 
physical space we begin to see the fronts evolving into highly convoluted sheets and the wave-wave interactions result 
in a redistribution energy from large to small scales [HJ],|I2- Of course as small-scale structures are being created, 
we see that the relatively inviscid behaviour seen for t < 1.5s ceases and dissipation of energy begins to increase. 



B. Enstrophy decay, invariant PDFs and spectral scaling 

Proceeding to the focus of this communication, from Fig. 1 we notice a marked change in the behaviour of the 
enstrophy as the fronts become severely distorted. Specifically the enstrophy, which grew during the initial front 
formation and subsequent development, now decays in a fairly monotonic manner. Examining the enstrophy in detail 
(see the log plot in Fig. 4), immediately after the maximum we notice that the primary signature of the decay is 
exponential along with a secondary small amplitude modulation. In fact, Fig. 5 shows the vorticity field well after the 
enstrophy attains it maximum value (see the figure captions for the exact times of the snapshots); as is expected via the 
non-conservation of vorticity from @, we do not see successive mergers resulting in large-scale structures, but rather 
the vorticity field continues to consist of distinct blobs separated by sharp ridges of concentrated enstrophy dissipation. 

Motivated by a somewhat similar scenario in the decaying passive scalar problem [2^ , [2!| , [3(| — where the decay 
(of the passive field) is purely exponential — we consider the normalized variable X = uj/Q 1 ^ 2 , where Q =< uj 2 > 
and < • > denotes spatial averaging. From J3J) and after performing a spatial average, the equation governing 
< X 2n > is 



il Jl < X 2n >= [7Q1 + agQ 2 ] < X 2n > -agQ 1 * < X 2n ~ 1 -— > -(2n - 1) 7 Q < X 2n ~ 2 (VX) 2 > (5) 
2n at ox 

where Q\ =< (Vcj) 2 > and Q2 =< lo^: >. In the passive scalar problem, X = 4>/Q)J 2 where <j> denotes the passive 

field and Q p =< 4> 2 >. In that case the purely exponential decay of < <j) 2 > (and higher moments) led to 9< * t > = 
and consequently PDF(A) attained an invariant profile |2^|,j33- However in the present case, the moments in (JjjJ 
inherit the secondary modulation from < uj 2 > and fluctuate about a mean. In fact, as all the moments have the 
same temporal fluctuations, the entire PDF is expected to attain an invariant shape but will exhibit small shifts in 
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magnitude. The extracted PDFs — see Fig. 6 — are plotted in three groups. The upper panel focusses on early 
times (t < 7s) and shows the approach to self-similarity. Interestingly, this approach is characterized by a gradual 
decrease in intermittency, i.e. the PDFs in the initial stages of evolution are extremely fat tailed — either stretched 
exponentials or power laws keeping in mind the difficulty in distinguishing between these two functions |3l||.|30j|. 

The middle panel of Fig. 6 shows the PDFs for the interval t £ [12, 28]s, i.e. the time during which the enstrophy 
experiences a modulated exponential decay. The self-similarity is evident, in fact the PDF's are now characterized 
by a small Gaussian core and an exponential tail. The 9 field during this interval is shown in Fig. 7 — notice that 
even though there exists a background layering, the frontal structures that ride on this background are oriented quite 
randomly. In other words, irrespective of its direction a one-dimensional cut of the snapshots in Fig. 7 will encounter 
a frontal jump. In effect the picture that emerges is, as the frontal structures become unstable (t ks 1.5s), energy that 
was trapped at larger scales during the very early stages of evolution begins to fill out the entire available range of 
scales. This "filling out" proceeds — from t rs 2 — 7s — until one achieves a state wherein the kinetic energy spectrum 
is close to a fc~ 5 / 3 scaling while the potential energy continues to scale as A; -1 (though at slightly larger scales) - 
this can be seen in Fig. 8 which shows the kinetic and potential energy spectra from t « 7s onwards. Of course as 
is seen in Fig. 8, the spectra gradually shift downwards as the total energy is decaying, the diffusive roll-off extends 
to larger scales and the scaling becomes less distinct (especially in the kinetic energy) as the fields become progres- 
sively smoother. In spite of this change in the nature of the flow, the enstrophy continues to decay in a modulated 
exponential manner and the self-similarity in the normalized vorticity field persits throughout the interval t £ [12, 28]s. 

Regarding the scaling of the energy spectra, the phenomenology of stably stratified turbulence suggests the 
presence of two regimes [2(J (see especially the discussion in Section 3 of Phillips H||), an inertial regime wherein 
the kinetic energ y s cales isotropically as fc -5 / 3 H^,|2^ and a bouyancy regime wherein the scaling is anisotropic 
and goes as k~ 3 |ll|,|2^. At *he present stage, Re » 1 coupled with Fr ~ 0(1) implies a turbulent regime that 
is weakly constrained by bouyancy. This, along with the observed directional independence of the spectra, suggests 
the applicability of inertial range phenomenology and thus supports the k~ 5 ^ 3 scaling of the kinetic energy. Given 
the weak bouyancy constraint, the scaling of the potential energy spectrum is dictated by the geometry of the 9— 
field. Specifically, as mentioned earlier even though the fronts become unstable they retain much of their identity, 
i.e. instead of being smoothed out, they evolve into highly distorted sheets (Fig. 7). In fact, these step-like features 
continue to dominate the potential energy spectrum resulting in the scaling [24[, p5|. 



C. Vertically sheared horizontal flows (a Pearson-Linden regime) 

Returning to the energetics of the flow, Fig. 1 indicates that at long times, i.e. t > 40s, the total energy reverts to 
an extremely slow decay. Further, in this stage almost all of the energy in the system is in the potential component 
(see also the third column of Fig. 2). The minute amount of kinetic energy indicates that a dissipative linear analysis 
would be appropriate — in fact, we are precisely in the last stages of decaying stratified turbulence as elucidated 
by Pearson & Linden p(l | (their analysis was more detailed with diffusivity ^ viscosity). Substituting a Fourier 
decomposition into the dissipative linearized form of 0, we have 

/ A _ 
ipik) exp{i(k x x + k z z)} exp{— "/K 2 t} exp{±i — \J (Xag)t} dk. (6) 
K 

The oscillatory nature of the exponential results in a dominance of the integral by modes with k x rj 0. Further, at 
long times a subset of these modes with the smallest rate of decay will remain. Hence we are left with a flow wherein 
w ~ 0, u — u(z) and 9 = 9(z), i.e. smooth vertically sheared horizontal flows with vertical structure restricted to the 
smallest wavenumbers (largest scales) as is seen in Fig. 9 ;40f. Regarding the vorticity field, from Fig. 4 we see that 
the nature of the enstrophy decay changes as we proceed into the Pearson-Linden regime — indeed, the set of PDF's 
in the lowermost panel of Fig. 6 correspond to these late times and as would be expected from a diffusive regime, the 
PDF gradually acquires more of a Gaussian form. 

Note that the form of © does not indicate how energy is transferred into the vertically sheared horizontal modes. 
Indeed a study of the forced 2D case [3^ suggests an inverse tranfer of energy into these quasi-horizontal modes, also 
Bartello ^ij^ demonstrates similar behaviour in the geostrophic energy during the process of geostrophic adjustment in 
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the 3D Boussinesq system. A detailed examination of this process is in progress — at this stage we can only conjecture 
that, much like behaviour of unstratified rotating fluids |34j|. at long times it is the near resonant interactions that 
feed energy into the slow modes (cr(fc) = 0) via interactions between fast waves (a(k) > 0) i.e. a slow-fast-fast transfer. 



IV. SUMMARY AND CONCLUSIONS 



We have studied the evolution of large-scale density (temperature) perturbations in the two-dimensional, stably 
stratified Boussinesq equations. The advantage of starting from a state of rest is the observation of various stages 
through which the system naturally evolves as governed by the Froude and Reynolds numbers. Starting from a 
smooth profile, wc immediately observe the formation of sharp fronts resulting in a frontally dominated k^ 1 potential 
energy spectrum. Further, the fronts spontaneously evolve into highly convoluted sheets accompanied by a spectral 
re-distribution of energy that culminates in the establishment of a fc~ 5 / 3 kinetic energy spectrum. Given the rapid 
decay of energy at this stage, the establishment of the aforementioned scaling is followed by a gradual downward shift 
in the spectra, the scaling becomes less distinct as the fields become progressively smoother and finally there emerges 
a large scale, slowly decaying, vertically sheared, almost horizontal mode wherein most of the energy is trapped in 
the potential component — i.e. the Pearson-Linden regime. 

With regards to the vorticity, the early stages of front formation and energy re-distribution are accompanied 
by a rapid increase in enstrophy. Then as the kinetic energy scaling is established we see the onset of an almost 
monotonic decay of enstrophy. In particular, immediately after the maximum the decay is primarily exponential 
with a secondary small amplitude modulation. Examining the normalized vorticity field, motivated by an analogous 
scenario in the decaying passive scalar problem, shows it to be characterized by an invariant exponential PDF. This 
almost exponential decay and the associated invariant PDF persists till we enter the Pearson-Linden regime. Finally, 
given the strong diffusive influence, deep into the Pearson-Linden regime the PDF relaxes towards a Gaussian profile. 

An interesting aspect of this problem is the inviscid limit, i.e. 7 — > 0. If the active scalar system maintains its 
regularity i.e. |V0|, |Vu|, |Vu>| < 00, then as 7 — > we expect ^ = 0. However, this does not imply 7 < (Vw) 2 >^ 0. 
In fact, given that the flow will not decay, we conjecture that the decay of enstrophy and the associated invariant 
PDFs will be established but the flow will never enter the Pearson-Linden regime. Further, the potential energy 
spectrum is expected to follow a k~ l scaling due to the presence of fronts, whereas the scaling of the kinetic energy 
is problematic due to a lack of dissipation. In fact, it is quite possible to end up with equipartition leading to a pile 
up of kinetic energy in the largest available wavenumbers [35j . On the other hand, if the system loses its regularity 
then H < even in the limit 7^0, and should behave as in the presently studied situation with fixed 7 > 0. It 
is worth noting that issues of a similar sort arise in the classical problem of geostrophic adjustment in the shallow 
water equations — as pointed out by Killworth j3|| (see also Kuo & Polvani |37j) — wherein the energy deficit of 
the geostrophic state has implications for wave-breaking and the generation of discontinuities from initially smooth 
data. In the present context it would be interesting to see how the difference in energy between the initial and final 
states scales with diffusivity. 

An obvious extension of the present work is the consideration of the fully three-dimensional (and also possibly 
rotating) problem — as mentioned in the introduction, the issue of balance and the spontaneous generation of 
imbalanced (or fast) waves is an active area of work Indeed, the qualitative similarity to certain aspects 

of decaying geostrophic adjustment in the rotating stratified three-dimensional Boussinesq system — elegantly 
elucidated by Bartello |3^ | — further motivates such an endevour. 
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List of Figures 

• Figure 1 : Upper Panel : Potential, kinetic and total energy with time. Lower Panel : Enstrophy as a function 
of time. The simulation is carried out at 750 x 750 resolution with 7 = 5 x 10~ 4 . The other parameters are 
specificcd in the main text. Note that the simulation was carried out till t ~ 175s, this plot only extends upto 
t w 60s so as not to be dominated by the Pearson-Linden regime. 

• Figure 2 : Same as Fig. (??) but the different stages are split up. Upper Panel : Potential, kinetic and total 
energy with time in the different stages. The first column shows the relatively inviscid evolution (till about 
t « 2s) followed by the generation of small scale structures. The second column starts with well defined spectral 
scaling of the KE and PE and signals the onset of enstrophy decay. The third column shows the slowly decaying 
Pearson-Linden regime where most of the energy is in the potential component. Lower Panel : Corresponding 
enstrophy plots. 

• Figure 3 : Snapshots of the temperature field, with reference to the Fig. (??) these are at t = 1.39, 2.34 seconds 
respectively. The emergence of the fronts perpendicular to the ambient stratification is clearly evident. 

• Figure 4 : log(< u 1 >) Vs. time that clearly shows a modulated exponential decay immediately after the 
enstrophy attains its maximum. The initial portion of the curve for very small times has been omitted for 
clarity. Note that at long times, as one enters the Pearson-Linden regime the rate of decay slows and its nature 
is quite different. 

• Figure 5 : Snapshots of the vorticity field during the modulated exponential decay of enstrophy. Note the field 
consists of blobs of vorticity separated by ridges of intense enstrophy dissipation. We encourage the reader to 
view these images in color as a grey scale printout masks the sharp features. 

• Figure 6 : PDFs of the normalized vorticity field. The upper panel shows the approach to a self-similar profile, 
note the decrease in intermittency with time. The two bunches of curves in the lower panel consists of profiles 
evenly spanning t <E [12, 28] sec and t 6 [126, 173] sec respectively. The upper bunch represents the interval 
during which the enstrophy decay is exponential in character (Fig. ??) whereas the lower bunch shows the 
PDF's in the Pearson-Linden regime. 

• Figure 7 : Snapshots of the 9— field corresponding to Fig. (??). Note that even though there exists a background 
layering, the fronts riding atop this layering are oriented in a fairly random manner. In other words, irrespective 
of direction a one-dimensional cut of the 9— field will encounter a frontal jump. 

• Figure 8 : Spectra of the kinetic and potential energy from ID cuts of the 9, u, w fields. The cuts are made in the 
vertical and horizontal directions and the scaling is identical in both the cases (these plots are the averages of the 
two cases) . The two bunches of curves have been shifted for clarity, the upper bunch are the KE spectra with the 
dot-dash line showing a fc~ 5 / 3 spectrum. Similarly the lower bunch are the PE spectra with the corresponding 
dot-dash line showing a spectrum. The kinetic energy attains k~ 5 ^ 3 scaling by about t « 7s. This is clearly 
seen in the compensated spectrum plotted at t = 6.7s, after which we see the the gradual downward shift of the 
spectra, the extension of the diffusive roll-off to larger scales and the scaling becomes progressively less distinct. 
Specifically, the plots are evenly spaced in the interval t 6 [6, 10]s. 

• Figure 9 : The establishment of a quasi-horizontal temperature field as we enter the final stages of decay — i.e. 
a Pearson-Linden regime. Referring to Fig. (??) we see that almost all the energy in the system is now in the 
potential component. 



